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ABSTRACT 

We investigate the total flux density, spectral, polarization, and Faraday rotation 
variability of HST-1 in the M87 jet during the outburst from 2003 to 2007 through 
multi-epoch VLA observations at 8, 15, and 22 GHz. Contrary to the general case for 
blazars, the flux densities of HST-1 rise earlier at lower frequencies from radio to X-ray, 
and the spectra are softening with the growth of outburst, indicating that the newly 
emerging subcomponents within HST-1 have relatively steep spectra. In particular, 
the intrinsic EVPA varies monotonically by ~ 90° at the 3 wavebands during the 
period, and all but the stationary subcomponent in the eastern end of HST-1 move 
downstream superluminally deviating divergently from the overall jet direction, with 
the motion of the outmost subcomponent bending from one side of the jet axis to 
another. These strongly argue for the presence of helical magnetic fields around HST- 
1, which is also supported by the fact that the subcomponents might be accelerated in 
this region. The fractional polarization is relatively low in the rising stage, and in the 
decaying stage the polarization levels are almost comparable at the 3 wavebands. In 
view of the quite large RM values, Faraday rotation is expected to occur dominantly 
external to HST-1 in the decaying stage, which is well supported by the presence of 
diffuse emission around HST-1, and consistent with the scenario that RM decrease 
gets slower with time. 

Key words: polarization — galaxies: active — galaxies: individual (M87) — galaxies: 
jets — radio continuum: general. 



INTRODUCTION 

Great progress has been made in understanding jet accelera- 
tion and collimation. Evidence is found that the optical elec- 
tric vector positional angle (EVPA) of a newly formed fea- 
ture rotates while it moves through a region upstream of the 
radio cores in BL Lac and PKS 1510-089, which is attributed 
to helical magnetic fields accelerating a nd focusing the jet in 
that zone IjMarscher et al.ll2008l . |2010| ). In these two cases, 
the accelerating and focusing process lasts until the feature 
encounters a region of turbulence, where the magnetic fields 
become chaotic. Whereas further attempts to investigate the 
process are still required, the jets in quite a few sources 
are found to be ac celerated downstream of the core (e.g., 
Ijorstad et al.l 120051 '). and helical magnetic fields are argued 
to be present within the jets in many sources by identify- 
ing the Faraday rotation measure (R M) gradients across the 
jet downstream from the core fe.g.. IZavala fc Tavlorll2005l : 
iGabuzda et al.l 120081 ). It is therefore likely that large-scale 
helical magnetic fields exist in jets, and play an important 
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role in jet acceleration and collimation downstrea m of the 
core as well, as suggested bv lLvutikov et al.l (|2005l ). 

M87 is one of the closest radio galaxies with a distance 
of 16. 7 Mpc, which is active from radio to 7-ray waveband 
fe.g.. lAcciari et aLl boOQ). HST-1 in M87 jet was defined to 
be a feature or component separated by ~ 0.86" from the 
central black hole, which is resolved into several 'subcom- 
ponents' at mas resolution with the most east one station- 
ary and some other s moving down the jet superluminally 
(|Cheung et al.ll2007l ). It seems that new subcomponents are 
ejected frequently from the stationary one. The flux density 
monitoring at multi-wavebands shows that a major outburst 
of HST-1 occurred from 2003 to 2007, of which the radia- 
tion was regarded to be synchrotron from radio to X-ray 
(IHarris et al.ll2003l. 12003). Along with a 7-ray flare occurred 
in the period (jAharonian et al.ll2006l ). HST-1 was argued 
to sh ow significant blazar-like properties fe.g.. [Harris et al.l 
[2003), of which the polarization nature has yet been rarely 
investigated. 

Here we present the multi-epoch full polarization results 
of HST-1 in M87, which are derived from VLA observations 
at 8, 15, and 22 GHz in the same A array. §2 describes 
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observations and data reduction, followed by the results and 
discussion in §3. Finally, a summary is given in §4. 



2 OBSERVATIONS AND DATA REDUCTION 

We collected the multi-frequency full polarization data from 
NRAO Data Archive q which were observed with the VLA 
A array including 8 epochs spanning the major outburst pe- 
riod from 2003 to 2007 to roughly keep consistency in reso- 
lution and sensitivity from epoch to epoch. The observations 
utilized two adjacent 50 MHz wide intermediate frequencies 
at all the three wavebands. 

Data calibration were performed in AIPS 
iBridel fc GreisenI 1 1994 ) with flux density and EVPA 
calibrated using 3C 286. Phases were initially calibrated 
through a nearby point source 1224-1-035. Instrumental 
polarization and EVPA corrections were obtained by 
using 'PCAL' and 'RLDIF', respectively. After the initial 
amplitude and phase calibration, th e data were then read 
into DIFMAP (shepherd et al.lll994l ') for imaging, with the 
resultant full polarization images loaded back into AIPS 
for model fitting using the task 'JMFIT'. Since HST-1 is 
not well separated from the core at 8 GHz , the whole 
core region containing HST-1 is taken to be fitted with 3 
elliptical components at 8 GHz, and 4 elliptical components 
at 15 and 22 GHz, and at all epochs the feature HST-1 is 
unambiguously identified. Note that since HST-1 is almost 
point-like at the VLA resolution, the effect of different beam 
size at the 3 frequencies on the results could be neglected. 
When model-fitting to the Stokes 'Q' and 'U', we fixed the 
size and position of each component derived from fits to the 
total intensity image to ensure that the resultant polarized 
intensity came from the same emission region (see Fig. [l]). 



3 RESULTS AND DISCUSSION 

The imaging results show that the jet extends to the north- 
west with a counter lobe structure clearly detected on the 
opposite side, and hence the detected emission region is as 
large as more than 50 arcseconds along the overall jet direc- 
tion at 8 GHz. In this Letter, we focus only on the feature 
HST-1, which changes significantly in various aspects dur- 
ing the major outburst. Figure [T] shows, just as an example, 
the central emission region of total intensity radio structure 
with fractional polarization (color) and EVPA vectors over- 
laid at epoch 2004 Dec. 31 at 8, 15, and 22 GHz. One can 
find that the total intensity radio structures appear similar 
at these frequencies with a little more extended structure 
detected at low frequencies. The fractional polarization and 
EVPA changes with distances to the core of M87. 

HST-1 is defined here to be the point-like emission 
feature excluding the extended structure downstream from 
HST-1, as is shown in Fig. [T] When the emission region up- 
stream of HST-1 is not taken into account, the radio struc- 
ture looks like a core-jet structure extending to northwest. 
The EVPA of HST-1 changes apparently with frequencies, 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Figure 1. Total intensity contours for the central emission re- 
gion of M87 with fractional polarization (color) and EVPA vec- 
tors overlaid at epoch 2004 Dec. 31. The contour levels start at 
1.5, 1.0 & 1.0 mjy, and increase by factors of 2, at 8 (X band), 15 
(U band), and 22 GHz (K band), respectively. The correspond- 
ing restoring beams are (0.26 X 0.24) arcsec'^ at P. A. = —7.6°, 
(0.16 X 0.14) arcsec^ at P.A. = 19.2°, and (0.10 X 0.09) arcsec^ at 
P. A. = 3.3° respectively, denoted as cross in the left of each image 
with its size indicating the beam size. The locations of HST-1 lie 
at places where a dashed line passes through at the 3 wavebands. 
The dashed-line box denotes the fitting emission region. 



which is due to F araday rotation f r om HST-1 itself an d its 
ambient medium (jChen et al.ll2010l : [iloman et al.ll2009 ll. All 
quantitative analysis with epoch and frequency, including 
the total flux density, polarization and RM of HST-1, are 
given in the following. 



3.1 Flux and Spectral Variability of HST-1 

The total flux density variations at the three wavebands are 
shown in Fig. [2^, of which the 15 GHz flux densities are in 
agreem ent with the results for the first several epochs avail- 
able in [Harris et al.l (|2009l ) . It is shown that the outburst 
grows relatively early at low frequency in the rising stage, 
which can also be f ound in a much w ider frequency range 
from radio to X-ray (JHarris et al.ll2009l ). This is contrary to 
the general case for blazars, where the flares at lower fre- 
quencies are gen erally delayed ownin g to synchrotron self- 
absorption fe.g:.. lMarscher et al.ll2008h . 

Since the flux density is relatively high at low frequency 
at all epochs (Fig. [TJi), the data is fitted over the 3 wave- 
bands in a power law spectrum at each epoch. One can find 
that the spectrum softens with epoch in the phase of rising 
flux, reaches the softest state in the peak stage, and then gets 
hardening in the decaying stage, as is shown in Fig.[2j3. Dur- 
ing the rising stage, at least two subcomponents are emerg- 
ing from the stationary subcomponent in the eastern end of 
HST-1 at times shown by arrows in Fig. [2^ IjCheung et al.l 
[20071), which implies that the newly emerging subcompo- 
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Figure 2. (a) Radio light curves of HST-1 at 8 (solid line), 15 
(dashed line) and 22 GHz (dotted line), the arrows indicate the 
times at which two subcomponents are emerging from the sta- 
tionary one within HST-1 on mas scale derived by Cheung et al. 
(2007); (b) Spectral variability of HST-1 with indices obtained 
through fits to the 3 waveband data. The errors to total flux den- 
sity are obtained in the course of model fitting, and the errors to 
spectral indices are derived in the course of fits. 



nent(s) has a relatively steep spectrum, and hence changes 
the hardness of the overall spectrum of HST-1 during the 
outburst assuming that the spectrum keeps constant for the 
stationary subcomponent. 

Note that the bright components of the radio jet appear 
to accelerate with distance from core to HST-1 with sublu- 
minal speeds, whereas several knots placed further out from 
the core (0.8-6.3 arc sec) move superluminally down the jet 
l|Stawarz et al.ll2006r ). In particular, multi-epoch optical and 
Very Long Baseline Array (VLBA) observations all showed 
that the unresolved stationary feature at the upstream edge 
of the HST-1 knot seems to emit vari ous subcomponents 
down the jet with relativistic speed (jBiretta et al.l 1 1999 : 
ICheung et al.ll2007l ) . That is to say, except for the stationary 
subcomponent which is special due to either standing shock 
or optical opacity there, all other subcomponents may well 
be accelerated, and electrons are energized in this region. 
On the other hand, a major flare of HST-1 occurred during 
the period of 2003 to 2007, accompanied by two subcompo- 
nents emerging from the stationary on e in the rising stage 
l|Cheung et al.ll2007l : [Harris et al.ll2009l) . Adiabatic compres- 
sion at least should not be the dominant mechanism for en- 
ergizing the radiating electrons, or else the spectr a would 
be hardening in the rising stage IjHarris et al.l 120031 ). which 



is contrary to t he observed r esults . It is more likely that, 
as suggested by [Harris et al.l (|2003l ). the increased intensity 
arises from the injection of new radiating particles that are 
energized in the region of HST-1 due to the fact that no 
signiflcant emission from radio to X-ray are detected im- 
mediately upstream of HST-1 even in the outburst period. 
Hence, there must be a way or mechanism, for example, a 
helical magnetic fleld, that plays a role in feature acceler- 
ation and particle energizing in the major outburst, which 
will be further discussed bellow. 



3.2 



Variations of Intrinsic EVPA and Fractional 
Polarization 



According to theoretical models (e.g., iMeier et al.l 1200 ll ). 
jets in AGN are propelled and coUimated by magnetic flelds 
twisted by differential rotation of the black holes accretion 
disk or inertial-frame- dragging ergosphere, and the primary 
observational indicator of magnetic coUimation requiring a 
helical magnetic fleld in the spine of the jet is the evolu- 
tion of the polarization, which manifests as EVPA rotation 
as a feature accelerat es and focuses along its spiral path 
JMarscher et al.ll2008l ). The EVPA rotation has been de- 
tected in the monitoring of EVPA to the so urces BL Lac and 
PKS 1510-089 (Marschc r et al.ll2OO8l , [2O10l ). Through muhi- 
epoch polarized VLA observations at 3 wavebands over 3 
years, we find a similar case with the feature HST-1 in M87 
jet. The intrinsic EVPA changes with time in a regular way, 
as is shown in the top panel of Fig. |3l In this map, the Fara- 
day rotation effect, which will give rise to EVPA change by 
Ax = RMA'^, has been eliminated at all the epochs. Here, 
RM is obtained through fits of EVPA to the wavelength 
square A^. Except for the epoch 2007 Aug. 10, at which the 
observational EVPAs of calibrator 3C 286, and hence M87 at 
22 GHz are not available, the resultant RM-corrected EVPA 
shows a roughly monotonous variation from ~ —90° to ~ 0°. 
Although the monotonous variation of EVPA may be also 
interpreted as a random walk of the resultant polarization 
vector direction as simulation cells with random magnetic 
field orientations enter and then exit the emission region, 
the p robability is supposed to be very low (jMarscher et al.l 
[2OI0I). It is more likely that the observed long-term EVPA 
va riation is dictated by th e magnetic geometry, as suggested 
bv lMarscher et al.l (2008). 

Another point is that as mentioned above, the emis- 
sion features may well be accelerated in HST-1 region, 
which is also supported by the fact that almo st all subcom- 
pone nts within HST-1 move superluminally (jCheung et al.l 
[20071) , and acceleration obviously occurs for HST-1 itself 
over the whole major outburst period (|Giovannini et al.l 
l2011f ). w hereas all componen ts upstream of HST-1 are sub- 
luminal (^Stawarz et al. 20061) ■ Based on the model proposed 
by [Marschcr ct al. ( 2008) , the emission features will move 
following a spiral path if a helical magnetic field exists, and 
different features may follow different subset of streamlines 
through the acceleration and coUimation zone. This scenario 
has been definitely detected by the multi-epoch VLBA ob- 
servations of M87 d uring the flaring period of HST-1 by 
ICheung et al.1 (|2007l ). The VLBA observations show that the 
resolved subcomponent a moves downstream at ~ 2.49c at 
P. A. ~ 295° initially, and then deflect to the direction of 
~ 289° with apparent speed of ~ 1.41c, differing to the over- 
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Figure 3. Bottom panel: Fractional polarization of HST-1 with 
time at 8 (bottom), 15 (midle) and 22 GHz (top); Top panel: 
RM corrected EVPA of HST-1 with epochs at 8 (bottom), 15 
(middle) and 22 GHz (top). The magnitude of error bars are 3 
times of uncertainty obtained through transfer from stokes T, 
'Q', 'U', and their errors. 



all jet direction of 290°. In addition, the other two subcom- 
ponents b and c move down the jet at P. A. ~ 279°, and the 
routes also appear not straight. It can be interpreted that 
different emission features may gain energy across different 
part of its cross-sectional area, and trace down the different 
subset of streamlines, and the changes in apparent speed 
ma y mainly attribut e to th e orientation change down the 
jet. iGiovannini et al.l (|201ir ) presented an e-EVN image at 
epoch 2010 Jan. 27, which shows a significant orientation dif- 
ference of the sub structures within HS T-1 to that at epochs 
by 2007 given in ICheung et all (|2007l '). If the helical mag- 
netic field around HST-1 do exists, each subcomponent may 
follow different subset of spiral streamlines down the jet, and 
with time, result in the significant orientation change of the 
substructures within it. All these effectively strengthens the 
argument that a helical magnetic field around HST-1 plays 
an important role in these phenomena. The approximate 100 
pc scale of the di sk of ionized gas observed in M87 by HST 
l|Ford et al.lll994l ) may indicate to some extent the presence 
of helical magnetic fields around HS T-1, which is also ab out 
100 pc away from the active centre (jStawarz et al.ll2006l l. 

Fractional polarization reflects the degree of magnetic 
field order, which may be affected by Faraday rotation 
and the medium opacity in the passage to the observer 
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Figure 4. Rotation measure of HST-1 with time. The top right 
inset is a A^ regression to the observed Faraday rotation at epoch 
2003 Aug. 24 with abscissa A^ in cm^ and ordinate EVPA in deg. 
The magnitude of error bars in the inset are 3 times of uncertainty. 



(jHoman et al.ll2009l ). and polarization level variability is ar- 
gued to arise from the s ource itself as well as its immedi- 
ate surrounding medium (|Homan et al.ll2009r ). We show the 
variations of the fractional polarization in the bottom panel 
of Fig. O which shows that the fractional polarization varies 
in a quite complex way. However, it is quite affirmative that 
the fractional polarization roughly minimizes in the peak- 
ing stage, and is relatively low in the rising stage, and high 
in the decaying stage. That is to say that the depolariza- 
tion effect due to Faraday rotation and medium opacity is 
relatively strong in the rising stage of the major outburst. 



3.3 Faraday Rotation of Emission from HST-1 

Faraday rot ation is often dete c ted from jet emission in radio- 
loud AGNs. lZavala fc Tavloj l|2004 l2005h suggested that it 
occurs external to the jet by excluding several other identi- 
ties such as broad and narrow line regions, and giving ev- 
idence of high fractional polarization to rule out internal 
Faraday rotation in 3C 273, while internal Faraday rota- 
tion gets supported in model-fitting to the full polarization 
spectra of 3 jet componen ts simultaneo usly in a consistent 
physical picture in 3C 279 (|Homan et al. 2009 ) . The rotation 
measure can be expressed as RM oc J N(s)B ■ ds, where N is 
the electron number density, B is the magnetic field, and s is 
the emission path to the observer. Obviously, RM may pro- 
vide constraints on number density and magnetic field in the 
passage, which is here obtained by fitting the observed EV- 
PAs over 3 wavebands with Faraday wavelength square law. 
The RM variation with time for HST-1 is shown in Fig.|4l in 
which the top right inset gives the fit of the observed EVPAs 
in deg to the squared wavelength in cm^, as an example, at 
epoch 2003 Aug. 24. One can find that the EVPAs obey the 
A^ Faraday rotation law well at the epoch. This means that 
although HST-1 contains complex substructures, its EVPA 
still obeys the Faraday rotation law well. 

Figure [4] shows that RM has the largest value in the 
beginning, and gradually decreases with fluctu ation present 
due to subflares occurred in the mean time IjHarris et al.l 
[2OO9I). This is consistent with the fact that fractional polar- 
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ization is relatively low in the early stage of the outburst, 
and as expected, the RM variability may well arise from the 
interior of HST-1. Faraday rotation may give rise to depo- 
lariza tion by a factor of sin(<E>)/cl? |Burn 1966; Homan et alj 
I2OO9I ). where $ = 2Ax, if Faraday rotation occurs internal 
to HST-1. According to the resultant RM values shown in 
Fig. [4] the dispersion effect is supposed to be quite signifi- 
cant at 8 GHz in comparison to that at the other two fre- 
quencies if Faraday rotation occurs internal to HST-1, for 
example, the polarization level at 8 GHz would be 4 times 
lower than that at 15 GHz even when RM took a relatively 
small value of 1000 rad m~ . Moreover, the relatively large 
beam size at low frequencies would make the fractional po- 
lariza tion even low er if it had an effect on the polarization 
level iJTavloii 119981 ). However, one can find in the bottom 
panels of Fig. [3] that after 2006, the polarization levels are 
almost comparable at the 3 frequencies. This implies that, 
although the RM variation is attributed to the plasma inter- 
nal to HST-1, quite a large part of Faraday rotation may oc- 
cur external to HST-1. Especially in the decaying stage, the 
Faraday rotation may occur dominantly external to HST-1, 
as is required by the comparable fractional polarization at 
the 3 wavebands. This is consistent with the fact that RM 
decrease gets slower with time, and more importantly, is 
supported by the presence of diffuse emission around HST-1 
m easured with VLA a t low frequencies, as is shown in Fig. [l] 
of [Harris et all (|2008l ). 



4 SUMMARY 

We reduced 8 epoch polarimetric observations of M87 with 
VLA A array at 8, 15, and 22 GHz from 2003 to 2007. 
Due to the s ignificant blazar-like nature shown in HST-1 
IjHarris et al.1 ^2008). the flux density, spectral, polarization, 
and RM variabilities are investigated in the period of a ma- 
jor outburst, during which a 7-ray flare is detected to be 
associated with it. 

Contrary to the general case for blazars, the outburst 
of HST-1 grows relatively early at low frequency, which is 
consistent with the spectral variability that the spectrum 
is softening in the rising stage. This implies that the newly 
emerging subcomponents have relatively steep spectra, and 
hence make the overall spectrum softening. In addition, 
these newly e merging subcompone nts from 2003 to 2004 are 
superluminal (jCheune: et al.ll2007^ ■ and hence expected to 
be accelerated in this region due to the fact that all compo- 
nents upstream of HST-1 are subluminal, and become faster 
with distance from the core. 

The RM-corrected EVPA of HST-1 shows a 
monotonous variation from ~ —90° to ~ 0° at the 3 
wavebands over a period of more than 3 years, and in 
this period the subcomponents within HST-1 move down 
along the routes deviating from the overall jet direction 
with the motion of one subcomponent bending from one 
side of the large scale jet axis to another, and another 
two subcomponents deflecting signiflcantly from the jet 
axis. These give quite strong evidence for the presence of 
a helical magnetic field around HST-1, which also gets 
support from the fact that the emission features may get 
accelerated in this region. The fractional polarization varies 
with epochs, and is relatively low in the rising stage in 



comparison to that in the decaying stage. This implies 
that the internal Faraday rotation and medium opacity 
play an important role to the variation of polarization 
level during the outburst. Moreover, in the decaying stage 
after 2006, the polarization levels are almost comparable at 
the 3 wavebands, indicating that Faraday rotation occurs 
dominantly external to HST-1 in this phase, or else it would 
be much lower at 8 GHz than measured. 
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